A model for projectile fragmentation 



00 



G. Chaudhuri^*, S. Malliki and S. Das Gupta^ 

^Theoretical Physics Division, Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, 
Kolkata 700064, India 

"^Physics Department, McGill University, Montreal, Canada H3A 2T8 



^ I E-mail: * gargi@vecc.gov. in 



Abstract. A model for projectile fragmentation is developed whose origin can be traced 
back to the Bevalac era. The model positions itself between the phenomenological EPAX 
parametrization and microscopic transport models like "Heavy Ion Phase Space Exploration" 
(HIPSE) model and antisymmetrised molecular dynamics (AMD) model. A very simple impact 
parameter dependence of input temperature is incorporated in the model which helps to 
analyze the more peripheral collisions. The model is applied to calculate the charge, isotopic 
distributions, average number of intermediate mass fragments and the average size of largest 
cluster at different Zbound of different projectile fragmentation reactions at different energies. 



1. Introduction 

Projectile fragmentation is an important technique for studying the reaction mechanisms in 
heavy ion cohisions at intermediate and high energies. In heavy ion cohisions, if the beam 
energy is high enough, the participant-spectator scenario can be envisaged. The participant 
zone is highly excited whereas the projectile like fragment (PLF) with rapidity close to that 
of the projectile rapidity and the target like fragment (TLF) with rapidity near zero are 
mildly excited. The PLF has been studied experimentally, this being one of the tools for 
production and identification of exotic nuclei. A projectile fragmentation model is developed 
[U [2] which involves concepts of heavy ion reaction plus the well known statistical model 
of multifragmentation (Canonical Thermo dynamical Model) and evaporation. This model is 
computationally much less intensive than heavy ion phase-space exploration (HIPSE) model [3] 
and antisymmetrized molecular dynamics (AMD) [3] which are based on transport calculation 
and less phenomenological than EPAX [5] which is based on the empirical parametrization of 
fragmentation cross sections . An impact parameter dependent temperature profile has been 
developed in order to better account for the results at different Zbound ranges and also to confront 
with data from different projectile fragmentation reactions at different energies. The model is 
in general applicable and implementable above 100 MeV/nucleon. 

The organization of the article is as follows. In Sec. 2 we describe the theoretical formulation of 
the model where as the impact parameter dependence of temperature is explained in Sec. 3. Sec. 4 
contains the results obtained from theoretical calculation and comparison with experimental data 
of different projectile fragmentation reactions with different projectile-target combinations and 
varying projectile energies. Finally summary and conclusions are presented in Sec. 5. 



2. Formulation of Model 

The model for projectile fragmentation reaction consists of three stages: (i) abrasion, (ii) 
multifragmentation and (iii) evaporation. In heavy ion collision, if the beam energy is high 
enough, then in the abrasion stage at a particular impact parameter three different regions 
are formed: (i) projectile spectator or projectile like fragment (PLF) moving in the lab with 
roughly the velocity of the beam, (ii) participant which suffer direct violent collisions and (iii) 
target spectator or target like fragment (TLF) which have low velocities in the lab. Here we 
are interested in the fragmentation of the PLF. The number of neutrons and protons in the 
projectile spectator at different impact parameters are determined from abrasion stage. Then 
the break up of each abraded projectile spectator is separately calculated by using canonical 
thermodynamical model (CTM) [6j. Finally, the decay of excited fragments are calculated by 
evaporation model [7] based on Weisskopf's formalism. We describe the details of the three 
different stages below. 

2.1. Abrasion 

In abrasion stage we assume the beam energy is high enough so that using straight-line geometry, 
PLF, TLF and participant can be classified. We then calculate the volume of the projectile that 
goes into the participant region (eqs. A. 4. 4 and A. 4. 5 of [8j). What remains in the PLF is 
V . This is a function of 6, the impact parameter. If the original volume of the projectile is 
Vq, the original number of neutrons is Nq and the original number of protons is Zq then the 
average number of neutrons in the PLF is < Ns{b) >= {V(b)/Vo)NQ and the average number of 
protons is < Zs{b) >= {V{b)/Vo)Zo. These will usually be non-integers. Since in any event only 
integral numbers for neutrons and protons can materialise in the PLF, we have to guess what 
is the distribution of Ng^Zg which produces these average values. For this purpose minimal 
distribution model is used. Let < Ns{b) >= N^^^{b) + a where a is less than 1. We can 
also define TVf "^(6) = iV™"(6) -Fl. We assume that Pat (6) is zero unless N,{b) is N'^''^{b) or 
jSi^ax^T^y The distribution is narrow. We then get P(iVJ""^'(^)) = a and P{N^^''{b)) = 1 - a. 
From < Zg > we can similarly define Pzs{b)- Together now we write Pns,Zs[^) = PNs{^)Pzs{b)- 
The abrasion cross-section when there are Ng neutrons and Zg protons in the PLF is labelled 
by (ya,N,,Zs'- 

(^a,Ns,Z. = J 2MbPN,,zAb) (1) 

where the suffix a denotes abrasion. The limits of integration in eq.(l) are bmin and bmax = 
Rtarget + Rprojectile- For bmin we have either (if the projectile is larger than the target) or 
Rtarget — Rprojectile (if the target is larger than the projectile, in this case at lower value of h 
there is no PLF left). Actually there is an extra parameter that needs to be specified. The 
complete labeling is (Ta,Ns,Zs,T if we assume that irrespective of the value of 5, the PLF has a 
temperature T. Here we have broadened this to the more general case where the temperature 
is dependent on the impact parameter b. In evaluating eq.(l) we replace integration by a sum. 
We divide the interval bmin 

to 

bmax into small segments of length A6. Let the mid-point of the 
i-th bin be < 6j > and the temperature for collision at < 5^ > be Tj. Then 

'ya,Ns,Zs =^(^a,Ns,Zs,Ti (2) 
i 

where 

cTa,7V„z„r, = 2t: <bi> AbPNs,Zsi< >) (3) 
PLF's with the same Ng, Zg but different Tj's are treated independently for further calculations. 



2. 2. Multifragmentation 

The abraded system of Ng neutrons and Zg protons created at impact parameter b will have an 
excitation which we characterize by a temperature T. The impact parameter dependence of the 
temperature profile is obtained from different projectile fragmentation reactions with different 
projectile target combinations. The details of it is described in the next section. The abraded 
system with Ng, Zg and a temperature T will break up into many composites and nucleons. We 
use the canonical thermodynamic model (CTM) [B] to calculate this break up. Assume that 
the system with Ng neutrons and Zg protons at temperature T, has expanded to a higher than 
normal volume and the partitioning into different composites can be calculated according to the 
rules of equilibrium statistical mechanics. In a canonical model, the partitioning is done such 
that all partitions have the correct Ng, Zg. The canonical partition function is given by 

"AT.Z 

Here the sum is over all possible channels of break-up (the number of such channels is enormous) 
which satisfy Ng = J2 ^ ^ nN,z and Zg = ^Z x un^z'-, ^n,z is the partition function of one 
composite with neutron number N and proton number Z respectively and n^^z is the number of 
this composite in the given channel. The one-body partition function lo^^z is a product of two 
parts: one arising from the translational motion of the composite and another from the intrinsic 
partition function of the composite: 

UJN,Z = ■^(27rmr)3/2^3/2 ^ ZN,ziint) (5) 

Here A = N -\- Z is the mass number of the composite and V is the volume available for 
translational motion; V will be less than Vf, the volume to which the system has expanded at 
break up (freeze-out volume) . We use V = Vf — Vq , where Vq is the normal volume of nucleus 
with Zg protons and Ng neutrons. In the projectile fragmentation model we have used a fairly 
typical value Vf = SVq. 

The probability of a given channel P{nN,z) = -P(?T'0.1) ^^1,0) '^i.i nj^j ) is given by 

P{nN,z) = -^ (6) 

The average number of composites with N neutrons and Z protons is seen easily from the above 
equation to be 

_ QNs-N,Zs-Z 

nN,z — ^N,z — (7) 

Qns,Zs 

There are two constraints: Ng = ^N x un^z and Zg = ^ Z x un^z- Substituting eq.(7) in 
these two constraint conditions, two recursion relations [9] can be obtained. Any one recursion 
relation can be used for calculating Qns,Zs- For example 

QNs,Zs = tt X! ^^n,zQns-n,Zs~z (8) 

N,Z 



We list now the properties of the composites used in this work. The proton and the neutron 
are fundamental building blocks thus zi^(){int) = zo^i{int) = 2 where 2 takes care of the spin 
degeneracy. For deuteron, triton, ^He and '^He we use ZN,z{i'nt) = {2sn,z + 1) exp(— /3£'7v,z(5'^)) 



where /? = l/T, i?7v,z(5?') is the ground state energy of the composite and {2sn,z + 1) is the 
experimental spin degeneracy of the ground state. Excited states for these very low mass nuclei 
are not included. For mass number ^4 = 5 and greater we use the liquid-drop formula. For nuclei 
in isolation, this reads {A = N + Z) 

ZN,ziint) = exp^[WoA-aiT)A'/'-a:^-CsyJ-^^^^ + ^] (9) 

The expression includes the volume energy [Wq = 15.8 MeV], the temperature dependent surface 
energy [a{T) = aoUT^ - T^)/iT^ + T^)}^/'^ with ctq = 18.0 MeV and = 18.0 MeV], the 
Coulomb energy with Wigner-Seitz approximation [TU] [a* = ac{l — (^o/^/)^^^} with Oc = 0.72 
MeV] and the symmetry energy {Cgym = 23.5 MeV). The term (eo = 16.0 MeV) represents 
contribution from excited states since the composites are at a non-zero temperature. 

We also have to state which nuclei are included in computing Qns,Zs (eq-(8)). For N, Z, (the 
neutron and the proton number) we include a ridge along the line of stability. The liquid-drop 
formula above also gives neutron and proton drip lines and the results shown here include all 
nuclei within the boundaries. 

We repeat the entire break up calculation for each projectile spectator created after abrasion 
stage with different temperatures at different impact parameters. Let, n^"'^"''^^ be the average 
number of fragment having A'^ neutron and Z proton created after the multifragmentation of 
a projectile spectator {Ns,Zs) at temperature Tj, then cross-section after multifragmentation 
stage can be expressed as 

(^m,N,Z,Ti = nN°,z"^'(^a,Ns,Z,,Ti (10) 

Ns,Zs 



2.3. Evaporation 

The excited fragments produced after multifragmentation decay to their stable ground states. 
Its can 7-decay to shed its energy but may also decay by light particle emission to lower mass 
nuclei. We include emissions of n,p,d,t,^lie and '^He. Particle decay widths are obtained using 
the Weisskopf's evaporation theory |TI]. Fission is also included as a de-excitation channel 
though for the nuclei of mass < 100 its role will be quite insignificant. According to Weisskopf's 
conventional evaporation theory, the partial decay width for emission of a light particle of type 
is given by 

r. = g^^o {E* - go - K) exp(2V^aMi^* - Eo - K) - 2 V^) (11) 

For the emission of giant dipole 7-quanta we take the formula given by Lynn [12] 

= r depR{E* - e)f{e) (12) 

pp[E*) Jo 

For the fission width we have used the simplified formula of Bohr- Wheeler given by 

Tf = ^exp{-Bf/Tp) (13) 

Details of the each term of the above three equations and implementation of the evaporation 
model can be found [7j; here we give the essentials necessary to follow the present work. Once 
the emission widths (F's) are known, it is required to establish the emission algorithm which 
decides whether a particle is being emitted from the compound nucleus. This is done [13] by first 



calculating the ratio x = r/rtot where Ttot = fi/^tot, ^tot = J2u ^"1/ and u = n,p, d, t, He^,a, 7 or 
fission and then performing Monte-Carlo sampling from a uniformly distributed set of random 
numbers. In the case that a particle is emitted, the type of the emitted particle is next decided 
by a Monte Carlo selection with the weights T^/Ttot (partial widths). The energy of the emitted 
particle is then obtained by another Monte Carlo sampling of its energy spectrum. The energy, 
mass and charge of the nucleus is adjusted after each emission and the entire procedure is 
repeated until the resulting products are unable to undergo further decay. This procedure is 
followed for each of the primary fragment produced at a fixed temperature and then repeated over 
a large ensemble and the observables are calculated from the ensemble averages. The number 
and type of particles emitted and the final decay product in each event is registered and are 
taken into account properly keeping in mind the overall charge and baryon number conservation. 
This is the third and final stage of the calculation. The same calculation is repeated for each 
set of fragments produced after multifragmentation at different temperatures. 

3. Temperature profile in projectile fragmentation 

Initially with the increase of the projectile beam energy, the temperature of the projectile spec- 
tator also increases. But above a certain energy of the projectile beam the temperature of the 
projectile spectator will not increase. This is known as limiting fragmentation. This projectile 
fragmentation model is valid in the limiting fragmentation region. The main reasons behind the 
excitation of projectile spectator are highly non-spherical shape and migration of some miclcons 
from participant to projectile spectator. Without a calculation at a more fundamental level it 
is not possible to calculate the excitation. At this stage we do not deal with excitation energy 
as such and characterize the system by a temperature T. We plan to work on this in future. 
Though the temperature is independent of incident energy of the projectile, but it depends upon 
the impact parameter. 

To get the impact parameter dependent temperature profile i.e. T = T{b) two types of 
parametrization comes to mind. The simplest case is that the temperature directly depends 
upon the impact parameter i.e. T(6) = Co + Cib + C2b'^ + ■■■■ In the another parametrization the 
temperature depends on the wound that the projectile suffers in the collision i.e. 1.0 — As{b)/Ao, 
so in this case T(6) = Do + Di{As{b)/Ao) + D2{As{b)/Ao)'^ + .... After calculating different ob- 
servables of projectile fragmentation by using temperature profiles and comparing the theoretical 
results with experimental data, it is observed that linear parameterizations are enough i.e. C2, 
C3... OT D2, D3... are negligible. Since T{b) = Dq + Di{As{b) / Aq) is physically more acceptable 
than T{b) = Cq + Cib, we finally choose this temperature profile and comparting with differ- 
ent sets experimental data at different projectile target combinations we select Dq = 7.5 and 
Di = —4.5. All theoretical calculations of this article is done by using this profile. 

For three different nuclear reactions ^^Ni on ^Be, ^^Ni on ""^^""^Ta and ""^^^Sn on ^""^^Sn, the vari- 
ation of the quantity Ag/Ao obtained after abrasion stage with normalized impact parameter 
{b — bmin)/ {bmax — bmin) is shown in Fig. l.a and Fig. l.b represents the freeze-out temperature 
profile of these three reactions calculated from the formula T{b) = 7.5 — 4.5(^5 (6)/Ao)- The 
specification that T(6) = Dq -\- Di{As{b)/Ao) has profound consequences. This means the tem- 
perature profile T{b/bmax) of ^^^Sn on ^^^Sn is very different from that of ^^Ni on ^Be. In the 
first case As{b)/Ao is nearly zero for b = &rmn=0 whereas in the latter case As{b)/AQ is ~ 0.6 for 
b = b.min=0- Even more remarkable feature is that the temperature profile of ^^Ni on ^Be is so 
different from the temperature profile of ®®Ni on ^^^Ta. In the latter case bmin = Rto — Rni and 
beyond bmin, -^s{b)/AQ grows from zero to 1 for bmax- This is very similar to the temperature 
profile of ^^^Sn on ^^^Sn. 
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Figure 1. Variation of (a) As/Aq and (b) temperature with normalized impact parameter 
(b - bmin)/ibmax - ^mm) for ^^Ni on ^Be (solid line), ^'^Ni on ^^^Ta (dotted line) and ^^^gn on 
^^^Sn (dashed line) reactions. 



In Fig. 2 the temperatures calculated from the model is plotted with Zf^ound {=Zs minus 
charges of all composites with charge Z = 1) and compared with experimentally measured 
temperatures (by Albergo formula [Ilj ) for two different projectile fragmentation reactions ^'^''Sn 
on ^^^Sn and ^^^Sn on ^^^Sn. The experiments are done by ALADIN collaboration in GSI at 
600A MeV il5j. 
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Figure 2. Comparison of theoretically used temperature profile calculated by the formula 
T{b) = 7.5— 4. 5(^5 (6) /ylo) (solid lines) with that deduced by Albergo formula from experimental 
data [15] (circles with error bars) for (a) ^''^Sn on ^^^Sn and (b) ^^^Sn on ^^^Sn. 



4. Results 

The projectile fragmentation model is used to calculate the basic observables of projectile 
fragmentation like the average number of intermediate mass fragments (Mjmf), the average 
size of the largest cluster and their variation with bound charge (Zbound), differential charge 
distribution, total charge distribution, isotopic distribution for different nuclear reactions at 
intermediate energies with different projectile target combinations. 



4.1. MjMF variation with Zbound 




Figure 3. Mean multiplicity of intermediate-mass fragments Mjmf, as a function of Zi^ound 
for (a) ^'^''Sn on ^^^Sn and (b) ^■^^Sn on ^^^Sn reaction obtained from projectile fragmentation 
model (solid lines). The experimental results are shown by the dashed lines. 



The variation of the average number of intermediate mass fragments Mjmf (3<Z<20) with 
Zbound for ^'^''Sn on ^^^Sn and ^^"^Sn on ^^^Sn reactions is shown in Fig. 3. The theoretical 
calculation reproduces the average trend of the experimental data very well. At small impact 
parameters, the size of the projectile spectator (also Zbound) is small and the temperature of the 
dissociating system is very high. Therefore the PLF will break into fragments of small charges 
(mainly Z = 1,2). Therefore the IMF production is less. But at mid-central collisions PLF's 
are larger in size and the temperature is smaller compared to the previous case, therefore larger 
number of IMF's are produced. With further increase of impact parameter, though the PLF 
size (also Zbound) increases, the temperature is low, hence breaking of dissociating system is very 
less (large fragment remains) and therefore IMF production is less. 

4-2. Differential charge distribution 

The differential charge distributions for different intervals of Zbound/ Zq are calculated by 
the projectile fragmentation model for ^^^Sn and ^'^^Sn on ^^^Sn reactions and compared 
with experimental data [15J. This is shown in Fig. 4. For the sake of clarity the 
distributions are normalized with different multiplicative factors. At peripheral collisions (i.e. 
0-8< Zbound /Zo<1.0) due to small temperature of PLF, it breaks into one large fragment and 
small number of light fragments, hence the charge distribution shows U type nature. But with 
the decrease of impact parameter the temperature increases, the PLF breaks into larger number 
of fragments and the charge distributions become steeper. The features of the data are nicely 
reproduced by the model. 
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Figure 4. Theoretical differential charge cross-section distribution (solid lines) for (a) ^'^''Sn 
on ^^^Sn and (b) ^^^Sn on ^^^Sn reaction compared with the experimental data (dashed lines). 



Size of largest cluster and its variation with Zi,ound 
Average size of the largest cluster produced at different Zf,ound values is calculated in the 
framework of projectile fragmentation model for ^^^Sn and ^^^Sn on ^^^Sn reactions. In Fig. 5 
the variation of Zmax/Zo {Z^ax is the average number of proton content in the largest cluster) 
with Zfjound/Zo obtained from theoretical calculation and experimental result are shown. Very 
nice agreement with experimental data is observed. 
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Figure 5. Zmax/Zo as a function of Zhound/Zo for (a) ^'^''Sn on -"^^^Sn (b) ^^^Sn on ^-"^^Sn 
reaction obtained from projectile fragmentation model (solid lines). The experimental results 
are shown by the dashed lines. 



4.4- Total charge distribution 

The total charge distributions of these experiments (^^Ni on ^Be and ^^Ni on ^^^Ta at 140 
MeV/nucleon done at MSU [16], ^^^Xe on '^^Al at 790 MeV/nucleon [H] and ^^^Xe on ^ospb 
at 1 GeV/nucleon |l8j both performed at GSI) are theoreticahy calculated from the projectile 
fragmentation model by using same temperature profile. This is shown Fig. 6. In theoretical 
calculation cross-section of all fragments ranging from light nucleon to original projectile are 
calculated separately, but in Fig. 6 the cross-section of the fragments for which experimental 
data is available is only shown. It is observed that, though the projectile beam energies in 
experiments are widely different, the same temperature profile can explain all the data pretty 
well. 
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Figure 6. Theoretical total charge distribution (solid lines) for (a) ^^Ni on ^Be, (b) ^^Ni 
on ^^^Ta, (c) ^^^Xe on ^^Al and (d) ^^^Xe on 208p|-, pg^ction compared with experimental data 
(dotted lines). 



4-5. Isotopic distribution 

For ^^Ni on ^Be reaction at 140 MeV/nucleon [T6], the isotopic distributions are theoretically 
calculated and compared with experimental values. This is shown in Fig. 7. Nice agreement 
between theoretical result and experimental data is obtained. 

5. Summary and Conclusion 

A model for projectile fragmentation is developed which is grounded on the traditional 
concepts of heavy-ion reaction (abrasion) as well as the model of multifragmentation (Canonical 
thermodynamical model) and secondary decay. This model is in general applicable and 
implementable in the limiting fragmentation region. An impact parameter dependent 
temperature profile is introduced in the model for projectile fragmentation which could 
successfully explain experimental data of different target projectile combinations of widely 
varying projectile energy. The observables which are calculated and compared to data included 
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Figure 7. Theoretical isotopic cross-section distribution (circles joined by dashed lines) for 
^^Ni on ^Be reaction compared with experimental data ^16j (squares with error bars). 



charge distribution, isotopic distribution intermediate mass fragment multiplicity and average 
size of largest cluster. 

While we have reasonable agreements with many data considered here, it is desirable to push 
the model improvements. The goal will be to find the size and excitation of initial projectile 
spectator from microscopic BUU calculations. We plan to work on this. 
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